Instant single-photon Fock state tomography 
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Heralded single photons are prepared at a rate of ~ 100 kHz via conditional measurements on polarization- 
nondcgencrate biphotons produced in a periodically poled KTP crystal. The single-photon Fock state is 
characterized using high frequency pulsed optical homodyne tomography with a fidelity of (57.6 ± 0.1)%. The 
state preparation and detection rates allowed us to perform on-the-fly alignment of the apparatus based on 
real-time analysis of the quadrature measurement statistics. © 2009 Optical Society of America 
OCIS codes: 270.5585, 270.5570, 270.5290 



The ability to produce single-photon states is of im- 
portance for optical quantum computation [1] . quantum 
cryptography [2], and quantum state engineering [3]. 
Many of these applications require the photons to be 
generated with a high fidelity in a pure, well-defined spa- 
tiotemporal mode. This requirement can be approached 
by heralded single photons prepared via conditional 
measurements on biphotons produced due to sponta- 
neous parametric down-conversion (SPDC) [4,5]. In this 
method, photon pairs produced in a nondegenerate op- 
tical parametric amplifier are split into two channels: 
trigger and signal. A single spatiotemporal mode is se- 
lected in the trigger channel and subjected to measure- 
ment with a single-photon counting module (SPCM). A 
detection event heralds preparation of the single pho- 
ton in the conjugate spatiotemporal mode of the signal 
channel [6,7]. 

A significant disadvantage of heralded photons is that 
they are not prepared on demand. The bimodal state of 
light produced due to each pump pulse (we specialize to 
the pulsed case) obeys thermal statistics: 

1*) « (1-72)(|0,0)+7|1,1)+7'|2,2) + ...), (1) 

where the numbers refer to Fock states in the signal 
and trigger channels and 7 is proportional to the pump 
field amplitude as well as the nonlinear susceptibility of 
the down-conversion medium. In the weak pump regime 
(7 <C 1), a trigger photon detection event indicates that 
the pulse in the signal channel contains a single photon. 
The probability of this event scales as oc 7^ . If multiple 
photons need to be generated from n heralded sources 
simultaneously, the probability scales as oc 7^", which 
results in unpractically long data acquisition times. A 
higher pump intensity increases the biphoton produc- 
tion rate, but leads, with non-ideal photon detectors. 



to a spurious multiphoton component in the signal, and 
thus to reduced fidelities. A compromise between these 
two regimes is facilitated by higher pump pulse repeti- 
tion frequencies. In this way, a reasonably high photon 
creation rate can be reached while maintaining a suffi- 
ciently low ratio between (n + 1)- and n- photon events. 
More frequent pump pulses however imply a reduced op- 
tical energy per pulse, which can be compensated by us- 
ing novel materials with high nonlinear susceptibilities. 
Additionally, the quantum state measurement optoelec- 
tronics must be sufficiently fast to distinguish between 
neighboring laser pulses. 

In this Letter, we report homodyne tomography [8] of 
the single-photon Fock state |1) generated under such 
intermediate conditions. The preparation rate is on the 
scale of 100 kHz with a pump pulse frequency of 76 MHz. 
This means, on one hand, that the two- and three-photon 
preparation rates are on the scale of 100 and 0.1 Hz, re- 
spectively, which is not prohibitively low. On the other 
hand, the fraction of multiphoton events heralded as sin- 
gle photons is on the order of 1% (assuming a 10% trig- 
ger photon detection efficiency), which is acceptable for 
most applications. A particularly attractive feature of 
our experiment is that the amount of data sufficient for 
reliable state reconstruction is acquired within a frac- 
tion of a second, permitting on-the-fiy optimization of 
the setup. 

Our setup can be compared with that described by 
Ourjoumtsev et at [9]. In that experiment, a cavity 
dumped Ti:Sapphire laser was employed to increase the 
pump intensity at the cost of decreasing the repetition 
frequency of the pump pulses to 800 kHz. Preparation 
rates of ~ 4500 Hz and ~ 15 Hz were reported for Fock 
states |1) and |2), respectively. 

Our apparatus is based on Ref. [6] and displayed in 
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Figure 1. (Color online) Experimental setup, fl: Spatial and 
spectral filters transmiting the pump light (blue) with A ~ 
395 nm. f2: Spectral filters transmitting states with A 790 
nm and AA = 0.3 nm. f3: lowpass filter blocking the pump 
light (Semrock BLP01-635R-25). 



Fig. [T] A mode-locked Ti:Sapphire laser (Coherent Mira 
900) emitted transform-limited pulses at ~ 790 nm with 
a repetition rate of 76 MHz and a pulse width of ~ 1.7 
ps. Most of the light was frequency doubled in a single 
pass through a 17 mm long LBO crystal, yielding typ- 
ically ~ 75 mW average second harmonic power after 
spatiotemporal filtering (fl). This field was coupled into 
a 2-mm long periodically poled potassium-titanyl phos- 
phate (PPKTP) crystal to pump type-II coUinear SPDC. 
The lengths of the crystals were chosen so that the group 
velocity dispersion would not lead to significant changes 
of the pulse width in both nonlinear transformations. 

The generated photon pair was separated into the sig- 
nal and trigger channels according to the photons' polar- 
izations using a A/2 plate and a polarizing beam splitter 
(PBS). A spatiotemporal mode for the trigger photons 
was selected with a 0.3-nm wide spectral filter (f2) cen- 
tered at the laser wavelength followed by a single mode 
fiber. The signal state was prepared in the conjugate 
mode when a SPCM (Perkin Elmer, SPCM-AQR-14-FC) 
registered the trigger photons. By temporarily blocking 
the pump field, we estimated that 2% of the trigger 
events occured due to scattered laser light, correspond- 
ing to a preparation efficiency of rjp ~ 0.98. Trigger event 
rates up to 140 kHz were routinely observed. 

A fraction of the original laser beam was split off be- 
fore frequency doubling to provide the local oscillator 
(LO) for balanced homodyne detection. The local os- 
cillator phase was allowed to vary randomly. The ho- 
modyne detector was aligned by coupling a 790 nm 
alignment beam into the trigger channel of paramet- 
ric down-conversion. Classical difference frequency sig- 
nal between the pump and alignment fields simulated 
the optical mode of the prepared signal photon [7, 10], 
and was matched to the LO by optimizing the interfer- 
ence visibility therewith. By implementing a zoom lens 
configuration in the LO path, a visibility of ~ 85% was 
reached, corresponding to a mode-matching efficiency of 
rjm ~ 73% [7]. The linear losses in the signal channel 
were mainly due to the lowpass filter (f3 in Fig. 1) and 



corresponded to an efficiency of rji ^ 96%. 

The homodyne detector circuit employed two Si-Pin 
photodiodes (Hamamatsu S5972) in a low-noise elec- 
tronic substraction circuit. These photodiodes possessed 
a quantum efficiency tjd ~ 85%, which is slightly lower 
than photodiodes typically used in high-speed homo- 
dyne detector circuitry in the past [8, 11-14]. However, 
their lower capacitance and dark current permitted us 
to significantly improve the bandwidth and signal-to- 
noise characteristics. The difference photocurrent was 
processed by an operational amplifier (OPA 847) in the 
transimpedance gain configuration. The detector had a 
—3 dB bandwidth at ^ 90 MHz, exceeding the repetition 
frequency of the LO pulses. The shot noise, measured at 
an average LO power of 25 mW, was ~ 14 dB higher 
than the electronic noise, corresponding to an electronic 
quantum efficiency of ijei ^ 93% [15]. 

The amplified difference photocurrent signal produced 
with the homodyne detector was digitized with an ac- 
quisition card (Agilent Acqiris DP211), integrated and 
displayed in real time. With each trigger event, 128- 
ns data segments were acquired. In this way, quadra- 
ture measurements were collected for the signal pulse as 
well as 8 neighboring pulses, whose state was approxi- 
mately vacuum, in order to calibrate the detector [15]. 
The acquisition rate, about 25,000 segments per sec- 
ond, was limited by the information transfer rate be- 
tween the acquisition card and the computer. The ac- 
quired quadrature data corresponds to the state p = 
(1 — 77) |0)(0| -f 77 11)(1|, where 77 is the overall quantum 
efficiency [6]. 

The variance of the acquired quadratures is then given 
by {Q^) = \ + V (using a scaling convention in which 
(Q^) = 5 for the vacuum state). The value of 77 ob- 
tained from the variance was calculated and displayed 
periodically on acquisition of a specific number of data 
segments, typically 4000 to 6000 (corresponding to a dis- 
play update rate of 4 to 5 times a second). The alignment 
of various parts of the setup could be improved using the 
displayed efficiency as a guide. Typically, starting from a 
preliminary alignment using classical fields, we were able 
to increase the efficiency by 5-10%. 

Upon optimization of the setup, a large data set of 
10® points was acquired and subjected to detailed anal- 
ysis [Fig. [2I^a.b)]. First, we corrected for residual corre- 
lations caused by a finite bandwidth of the homodyne 
detector, which lead to an additional quantum efficiency 
improvement of 0.5% on average. The maximum likeli- 
hood method [8,16,17] was then used to reconstruct the 
diagonal density matrix elements shown in the left col- 
umn of Table 1 and in Fig. [D^c). Contributions of higher 
number states were negligible. 

The vacuum state data yielded the density matrix el- 
ements displayed in the right column of Table 1 with 
negligible errors [Fig. [Hd)]. The presence of the two- 
photon component in the vacuum state estimation is due 
to the persistent non-Gaussian "hump" around Q = Q in 
the marginal distribution [Fig. [^b)] which we believe to 
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Figure 2. (Color online) Experimental results. Top row: his- 
togram of the collected signal state (a) and vacuum state (b) 
quadratures. The black lines are scaled marginal distributions 
based on the reconstructed density matrices. Bottom row: Re- 
constructed diagonal density matrix elements for the signal 
(c) and vacuum (d) states. The insets show cross-sections of 
the corresponding Wigner functions. 



Table 1. Quantum state reconstruction results. The uncer- 
tainties shown are standard deviations of the state estimation 
based on the Fisher information matrix [18]. 





Fock 


vacuum 


Poo 

P22 


0.4138 ±0.0006 
0.5758 ±0.0011 
0.0104 ±0.0006 


0.9987 


0.0013 



arise due to a finite digitizer resolution. 

The noticeable value for p2,2 = 0.0104 ± 0.0006 in the 
Fock state density matrix is likely to be caused by mul- 
tipair creation events in parametric down-conversion. In 
order to evaluate the suitability of our setup for mul- 
tiphoton quantum technology experiments, let us esti- 
mate 7. Assuming the quantum efficiency of the SPCM, 
including linear losses in the trigger channel, to equal 
r]t <^ 1, the state heralded by the photon detection 
event, neglecting dark counts, is given, according to 
Eq. O, hyph a mi^ |1)(1| + 2r/t7'* |2) (2| + 0(76). If this 
state is subjected to an optical loss 1 — 1]. it becomes 
p oc ??t[(l - 77)72 + 2(1 - 77)274] |0)(0| + 7747772 + 477(1 - 
77)74] |1)(1| + 27/47/274 |2}(2| + 0(7^). Normalizing and 
comparing this result with the density matrix obtained, 
we find 77 = 0.5787 ± 0.0011 and 72 « 0.0160 ± 0.0010. 

The probability of the heralding event is pt ~ ^7*7^- 
In this particular data set, the trigger event frequency 
was 88 kHz, which corresponds to pt = 1.15 • 10~^, from 



which we find 774 « 7%, which is consistent with direct 
estimate of the trigger channel transmission. 

In summary, we have reported homodyne reconstruc- 
tion of the single-photon Fock state produced using con- 
ditional measurements on photon pairs generated via 
type-II spontanteous parametric down-conversion in a 
periodically poled KTP crystal with a detection rate on 
the scale of 100 kHz. The single-photon reconstruction fi- 
delity is comparable to that reported previously [6,9,13]. 
The combination of the preparation rate as well as fast 
measurement and analysis permitted virtually instanta- 
neous tomographic characterization, so optimization of 
the experimental setup in real time was possible. The 
achieved parameters of the experimental setup permit 
implementation of various new experiments in quantum- 
optical engineering and information processing that re- 
quire simultaneous operation of two and three single- 
photon sources. 
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